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Abstract 



In reaction rate theory, in input-output type models, and in reaction-diffusion problems, when the 
total derivatives are replaced by fractional derivatives the solutions are obtained in terms of Mittag-Leffler 
functions and their generalizations. When fractional calculus enters into the picture, the solutions of these 
problems, usually available in terms of hypergeometric functions, G and H-functions, switch to Mittag-Leffler 
functions and their generalizations into Wright functions. In this paper, connections are established among 
generalized Mittag-Leffler functions, pathway model, Tsallis statistics, superstatisitcs, and power law, and 
among the corresponding entropic measures. 
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1. Introduction 

Fundamental laws of physics are written as equations for the time evolution of a quantity x(t), 

dj f = Mt), (i.i) 

where if A is limited to a linear operator we have Maxwell's equation or Schroedinger equation, or it could 
be Newton's law of motion or Einstein's equations for geodesies if A may also be a nonlinear operator. When 
A is linear then the mathematical solution is 

x(t) = x a c- At (1.2) 

where x$ is the initial value at t = 0. In reaction rate theory, if the number density at time t of the -i-th 
particle is Ni(t) and if the number of particles produced or the production rate is proportional to Ni(t) then 
the reaction equation is 

^ = hN(t), h > 

deleting i for convenience. If the decay rate is also proportional to N(t) then the corresponding equation is 

^ = -k 2 N(t), k 2 > 0. 
Then the residual effect in a production-destruction mechanism is of the form 



dN(t) 
dt 



= -cN(t), c> N(t) -N = -c N(t)dt (1.3) 



if the destruction rate dominates so that the input-output model is a decaying model. Such input-output 
models abound in various disciplines. If the total integral or the total derivative in (1.3) is replaced by a 
fractional integral then we have 

N(t) -N = -c v Dt U N(t), (1.4) 
where qD^ v is the Riemann-Liouville fractional integral operator defined by 

oDrm = -f- At - ur-\f(u)du (1.5) 

r W Jo 

for 3i(i/) > where Sft(-) denotes the real part of (•), c is replaced by c u for convenience, and then the solution 
of (1.4) goes into the category of a Mittag-Leffler function ([3]), namely 

"(«) = witSt = "° ^ [ - (c * rl (L6) 

fe=0 v ; 

where S^(-) is the Mittag-Leffler function. The generalized Mittag-Leffler function is defined as 

El,( Z ) = ± w ^ y ^a)>0, m >0, (1.7) 
and when T(7) is defined, it has the following Mcllin-Barnes representation 



KA Z ) = rw^Mffi + rt) (L8) 

- J- J- f ^ »l!l(-,)-Ma (1 9) 

- r( 7 ) 2™ y c _ ioo r(/3- as ) ( zj (L9) 
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for < c < SR(7),5R(7) > 0,i 
following: 



-1. Some special cases of the generalized Mittag-Leffler function are the 



oo 



k=0 



T(/3 + ka) 



E, 



k=0 



(1 + ka) 



(1.10) 
(1.11) 



and when a = 1 we have 



E a (z) = E 1 (z) = e z . (1.12) 

Thus the Mittag-Leffler function can be looked upon as an extension of the exponential function. The 
Mcllin-Barnes representation in (1.9) is a special case of the Mellin-Barnes representation of the Wright's 
function ([19], [20]) which is defined as 



{b j ,0 j ),j=l,...,q_ 



oo 



2niJ c 



c+ioo 



T(s) 



(-z)- s ds 



(1.13) 
(1.14) 



where < c < mmi<j< p 3? (^^j with cij,j = l,...,p and bj,j = l,...,q being complex quantities and 

ctj > 0, j = 1, ...,p and (3j > 0,j — 1, ...,q being real quantities. Observe from (1.14) that Wright's function 
is a special case of the H-function ([4], [10], [11]) and the H-function is defined as the following Mellin-Barnes 
integral 



H p,q \ Z ) = H p,q 



(oi,ai),...,(op,a p ) 
(6i,/3i),...,(6„/3,) , 



i pc+ioc 

— / <Ks) z -6s 

AH I Jc—ioo 



where 



<t>{8) = 



{ut=i r fe + m} {uu r (! - a i - 
{iTU+i r(i - 6^ - m} {i\-= n+ i r (^- + 



(1.15) 



(1.16) 



for maxi<j< m 3 ^ < c < mini<j<„ J?( ' 1 Q ^ a: '' > where a,j,j = 1, and &j, j = 1, q are complex quantities, 
Qj > 0, j = 1, and f3j > 0,j = 1, q are real quantities. Existence conditions and various contours may 
be seen form books on H- functions, for example, ([4], [10], [11]). 



1.1. Extension of the reaction rate model and Mittag-Leffler function 

The reaction rate model in (1.3) can be extended in various directions. For example, if N is replaced 
by N f(t) where f(t) is a general integrable function on the finite interval [0,6] then it is easy to see that 
for the solution of the equation 



N(t) - N f{t) = -c v Dt"N(t) 
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(1.17) 



there holds the formula 



N(t) = cW ( 



/ H l',2 

Jo 







c v (t~x) v 









f(x)dx. 



(1.18) 



Some special cases are the following: Let v > 0, p > 0, c > 0. Then for the solution of the fractional equation 

N(t) - Not"- 1 = -c v oD~ v N(t) 



there holds the formula 

N(t) = NoTW^E^i-icty). 
Let c > 0, v > 0, p > 0. Then for the solution of the fractional equation 

N(t) - Not^EZJ-W] = -c v D^N(t) 

there holds the formula 

N(t) = N Q t»- 1 EJ+ 1 [-(cty}. 
1.2. Fractional partial differential equations and Mittag-Leffler functions 



(1.19) 
(1.20) 
(1.21) 
(1.22) 



In a series of papers, see for example, ([14], [15], [16]) it is illustrated that the solutions of certain fractional 
partial differential equations, resulting from fractional diffusion problems, are available in terms of Mittag- 
Leffler functions. For example, consider the equation 



a D v t N(x,t) 



(1.23) 



with initial conditions 



Q D»- k N(x,t)\ t=Q =0,k= l,...,n 

where n — [SR(i^)] + 1, c u is the diffusion constant, 5(x) is the Dirac's delta function and [Jt(f)] is the integer 
part of 5i(z/). The solution of (1.23), by taking Laplace transform with respect to t nd Fourier transform 
with respect to x and then inverting, can be shown to be of the form 



N(x,t) = 



r/2,0 



(A-KC nu t u )* 

which in special cases reduce to Mittag-Leffler functions. 



r ixi 2 


(!-«.") " 




(0,1), 



(1.24) 



This paper is organized as follows: Section 2 gives the classical special function technique of getting rid 
of upper or lower parameters from a general hypergeometric series. Section 3 establishes the pathways of 
going from Mittag-Leffler function to pathway models to Tsallis statistics and superstatistics through the 
parameter elimination technique. Section 4 gives representations of Mittag-Leffler functions and pathway 
model in terms of H-functions and then gives a pathway to go from a Mittag-LcfHer function to the pathway 
model through H-function. 

2. A classical special function technique 

The age-old technique of getting rid off a numerator or denominator parameter from a general hyper- 
geometric function in the classical theory of special functions, is the following: For convenience, we will 
illustrate it on a confluent hypergeometric series. 



i F 1 (a;b;z) 



(a) 



a(a + l)...(a + m — 1), a ^ 0, (a) = 1. 



(2.1) 
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Observe that 



(a)k_a(a + l) (a + ra-1) 

a k a a a 

1 9 in — 1 

= 1(1 + -)(1 + -)...(1 + ) ^lasa^oo (2.2) 

a a a 

for all finite k. Similarly 

b k 

— > 1 when b — > oo. (2.3) 

Therefore 

lim l F 1 {a;b;-) = F 1 { ;b;z) (2.4) 

a^oo a 

which is a Bessel function. Thus we can go from a confluent hypergeometric function to a Bessel function 
through this process. Similarly 

lim 1 F 1 (a;b;bz) = 1 F (a; : z) = (1 - z)~ a , \z\ < 1. (2.5) 
Thus, from a confluent hypergeometric function we can go to a binomial function. Further, 



and 



lim Fi( ;6;6z) = F ( ; ; z) - e z (2.6) 
lim !F (a; ;-) = o*o( ; ;^) = e z . (2.7) 

a^oo a 



Thus, we can go from a Bessel function as well as from a binomial function to an exponential function. These 
two results can be stated in a slightly different form as follows: 

lim 0*i ( *—) = e~ z (2.8) 

g->l q — 1 (/ — 1 

and 

l imiFo (_L ; ;-( ( ,-l)z)=lim[l + ( ? -l) z ]-FT= c - (2.9) 
g-»l g — 1 g->l 

Equation (2.9) is the starting point of Tsallis statistics, non-extensive statistical mechanics and g-calculus. 
The left side in (2.9) is the g-exponential function of Tsallis, namely, 

[l + {q-l)z]-^ = exp q (-z) 

= c~ z when q -» 1. (2.10) 

We consider a more general form of (2.9) given by 

lim c\xPiF (-^— : ; -a(g - l)|a;| 4 ) = lim c|x| 7 [l + a(g - l)|a;|*]~^ 
9-»l q — 1 

= c|a;| 7 e- a7 'l x l 5 ,a > 0,j? > 0,-oo < .t < oo. (2.11) 
Thus, we obtain the pathway model of ([5]) for the real scalar case, namely, 

f(x) = c\x\~<[l + a(q-l)\x\ s }-^,a> 0,7? > (2.12) 
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where c is the normalizing constant. If r\ = 1, 7 = 0, o = 1, S = 1, x > then (2.12) gives Tsallis statistics 
([17, 18]). Hundreds of papers are published on Tsallis statistics showing a wide range of applications of the 
function. For q > 1, a = 1, 7/ = 1, x > in (2.12) we get superstatistics ([1],[2]). Dozens of articles are written 
on superstatistics showing the variety of practical situations where the concept of superstatistics conies in. 

It may be observed that (2.7), which is the binomial form or iFq going to exponential, is exploited to 
produce the pathway model, Tsallis statistics and superstatistics whereas (2.6), which is the Bessel function 
form going to exponential is not yet exploited. This should also produce a rich variety of applicable functions. 

3. Connection of Mittag-Leffler function to the pathway model 

In order to see the connection, let us recall the Mcllin-Barncs representation of the generalized Mittag- 
Leffler function. 

s m 1 r+- r( 7 - s )r( s ) s 

for 3^(7) > 0,2 = \f—l, $t(/3) > 0. Let us examine the situation when |/3| — > 00. We will state two basic 
results as lemmas. 

Lemma 3.1 When — > 00 and a > is finite then 

,i im T [-K*/?*)V = H»*V- (3-2) 

|/3|->oo 1 (p — as) 

Proof. The Stirling's formula is given by 

T(z + a) w V2^z z+a -ie- z , for |z| -» 00 (3.3) 
and a is a bounded quantity. Now, applying Stirling's formula we have 



r(/3-as) L w ;J y/2H0P-7- a8 e-P 

= I3 as {-bz 5 )- s l3- as = {-bz 5 )- s . (3.4) 



Lemma 3.2. For 3?( 7 ) > 0, 5R(/3) > 7 



1 1 rc+ioo 

lim r(f})EZj(-bF*z s ) = W ~^ r( 7 - s)r( S )[^]- s d S (3.5) 



|/3|-oo 



= ^ <^(-&/) fe = [1 + fo^H /or |6z 5 | < 1. (3.6) 

fc=0 

Proof. After applying Lemma 3.1 we can evaluate the contour integral in (3.5) by using the residue theorem 
at the poles of T{s) to obtain the series form in (3.6). For \bz s \ > 1 we obtain the series form as analytic 
continuation or by evaluating the integral in (3.5) as the sum of the residues at the poles of r( 7 — s). This 
will be equal to the following: 

lim T(l3)El (-bf3 a z s ) = (bz s )^[l + (bz*)- 1 }- 1 (3-7) 

for \bz s \ > 1 which will be of the same form as in (3.6). Replacing b by a(q — 1) and 7 by we have the 
following: 
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lim c\zpr(0)lFj i l°- 1 \-a(q-l)0 a \z\ s ) 

= c\zp[l + a(q-l)\z\ s }-^. (3.8) 

Observe that (3.8) is nothing but the pathway model of ([5]) for a > 0, S > 0, r\ > with c being the 
normalizing constant, for both q > 1 and q < 1. Then when q — > 1 (3.8) will reduce to the exponential form. 
That is, 

lim lim c|z|T(/3)^^ 1) [-a(9-l)/3 a kl' 5 ] = c|z|T e - a,)|z|i . (3.9) 

Note from (3.8) that when 9 > 1 the functional form in (3.8) remains the same whether S > or S < with 
—00 < z < 00. When (7 < 1 the support in (3.8) will be different for 5 > and <5 < if f(x) is to remain as 
a statistical density. 

For 7 = 0, a = 1, 5 = 1, 77 = 1, 2 > we have Tsallis statistics coming from (3.8) for both the cases q > 1 
and g < 1. For q > 1, a = 1, 6 = 1, 77 = 1, 2 > we have superstatistics coming from (3.8). Observe that 
the limiting process from (3.5) to (3.6) holds for 7 = 1 or 7 = 1 and a = 1 also. Thus the special cases of 
Mittag-Leffler function are also covered. 

Thus, through a pathway is created to go from a generalized Mittag-Leffler function to Mathai's 
pathway model and then to Tsallis statistics and superstatistics. If is real then as becomes larger and 
larger then 

c\zm0)E^- 1 \-a{q-l)0 a \z\ s ) 

goes closer and closer to the pathway model in (3.8). In other words, a pathway is created through to go 
from a Mittag-Leffler function to the pathway models to Tsallis statistics and superstatistics. Thus for large 
real value of or for large value of \0\, 

c\zyT(0)E 11 J^- 1) [-a{q - l)0 a \z\ 5 ] w C |zp[l + a(q - l)\z\ s }~^ . (3.10) 

4. Connections through the H-function 

Recalling the generalized Mittag-Leffler function from (3.1) and representing it in terms of a H-function 
we have the following: 



T{0)EZ^z0$Y = 



r(/3) 1 

r( 7 ) 2m 

r(7) 1>a 



c+»oo r ( s ) r ( 7 _ g) 

-too r(/3 - as) 

(1-7,1) 

-(z0T) s 

(0,1), (1-/3, a) 



[-(z/3f) 5 r s d S 



(4.1) 
(4.2) 



From the limiting process discussed in (3.2) we have the following result: 
Lemma 4.1 For 5ft(/3) > 0, 3?( 7 ) > 0, 



lim ——H, 9 

i/jHoor( 7 ) M 



^ J I (0,1), (1-/3, a) 



" IXt)^ 1 
1 1 /■ 

~ T(j)2mJ c _ 

= n - z'l 



.51(1-7,1) 



1(0,1) 



C+ZOO 



r( s )r( 7 - s )(-z 5 )- s ds 



(4.3) 

(4.4) 
(4.5) 
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Therefore we have the following theorem: 



Theorem 4.1 For SR(/3) > 0, ^(7) > 0, x > 0, a > 0, q > 1, c > 

lim c F(/3) . a;-y^/(g-i)r_ ( g _ 



|/3|—oo p 



(*)' 



ca; 7 [l + a(g- 



(4.6) 



Observe that the right side in (4.6) is the pathway model ([5]) for x > from where one has Tsallis statistics, 
superstatistics and power law where the constant c can act as the normalizing constant to create a statistical 
density in the right side of (4.6). 

We will write the right side in (4.6) as a H-function and then establish a connection between Mittag- 
Leffler function and the pathway model through the H-function. To this end, the practical procedure is 
to look at the Mellin transform of the right side of (4.6) and then write it as a H-function by taking the 
inverse Mellin transform. The Mellin transform of the right side of (4.6), denoted by Mf(s), is given by the 
following: 

/•OO 

M f {s)= I cx'> +s - 1 [l + a(q-l)x s }-^dx,a>0,q> 1, 77 > 0, 5 > 



(4.7) 



for 3? (-^p^ > 0, 3? i^pi — ~^T^j > 0- Note that if c is the normalizing constant for the density f(x) then 
by putting s = 1 the right side of (4.7) must be 1. Therefore, 



M f (8) = 



K«-i)]* 



T] 7 + S 



9-1 



rmr(i?i-f) V * 

x [o(«- I)] - *. 

Hence the right side of (4.6) is available as the inverse Mellin transform of (4.8). That is, 



(4.9) 



where 



f(x) = cx^[l + a(q- l)x 5 ]-^, 

for a > 0, 7] > 0, q > 1, 8 > 0, x > 
[a(g-l)]* 



77 7 + s 



9-1 



x j_ r +io ° T (i+jL) 

2™ 7 c „ loo V s J 
r(^)r(^-^i) 

= J™ «=-^ T ^J'- 1 »[- I ( B ( ! - !))*/?* ]« 

V 9-lJ 



[z(a(g- l))*]" s ds 



-if 



1.1 



x(a(«- 1))« 











(4.10) 



(4.11) 
(4.12) 

(4.13) 
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_ S[a( q -1)]^T(^) 

C_ r(2±i)r(^-2±i) 

for q > 1, S > 0, rj > 0, a > 0, R( 7 + 1) > 0, - > 0. 

Remark 4.1. From (4.13) and (4.10) it can be noted that as the parameter becomes larger and larger 
the Mittag-Leffler function in (4.13) comes closer and closer to the pathway model and eventually in the 
limiting situation both become identical. In a physical situation, if (4.10) represents the stable situation 
then the unstable neighborhoods are given by (4.13). When q — > 1 then (4.10) goes to the exponential form. 
Thus if the exponential form or Maxwell-Boltzmann situation is the stable situation then the pathway model 
of (4.1) itself models the unstable neighborhoods. This unstable neighborhood is farther extended by the 
Mittag-Leffler form in (4.13). 

Remark 4.2. Observe that the functional form in (4.10) remains the same whether q > 1 or q < 1. But 
for q < 1 or when q — > 1 the normalizing constant c will be different. 

Remark 4.3. When dealing with problems such as reaction-diffusion situations or a general input-output 
model one goes to fractional differential equations to get a better picture of the solution. Then we usually 
end up in Mittag-Leffler functions and their generalizations into Wright's function. Comparison of (4.13) 
and (4.10) reveals that as gets larger and larger the effect of fractional derivative becomes less and less and 
finally when \0\ — > oo the effect of taking fractional derivatives, instead of total derivatives, gets nullified. 
Physical interpretation of this pathway of going from (4.13) to (4.10) and then to exponential can hopefully 
produce new physics or explanations to currently unknown phenomena. 

5. Mittag-Leffler to Levy Distribution 

The generalized Mittag-Leffler density given by 

a/3-1 oo ro\ (- r a ) k 

f w - — £ { -f stU I an • * * < * > » > 

has the Laplace transform 



(4.14) 



(5.1) 



L f (t) = [l + 6t a }-e. (5.2) 
If 6 is replaced by 5(q — 1) and by 0/(q— l),q> 1 and if we consider q approaching to 1 then we have 

lim L f (t) = lim[l + %- l)t Q ]-^r = e' 801 " . (5.3) 

But this is the Laplace transform of a constant multiple of a positive Levy variable with parameter a, < a < 
1, with the multiplicative constant being (S0)~ , and thus the limiting form of a Mittag-Leffler distribution 
is a Levy distribution. Writing f(x) as a Mellin-Barncs integral we have the following: 



'<*> -mkhSZ W^T^ 0<c< " (5 - 4) 

1 1 j-ci+ioo — i + -)r(l + — — — ) X 

= fC5)2^i Cl _ ioo a S iT(2- s ) " a fe rMs ' (5 ' 5) 

Hence if f(x) is a density then we can take the kernel in the Mellin-Barnes integral as the (s— l)-th moment 
E^x 3 - 1 ). Thus 



Then s = 1 should give 1. The right hand side gives 1 and hence f(x) in (5.1) is a density function. This is 
called the generalized Mittag-Leffler density. 
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